Transient magnetoconductivity of photo excited electrons 

O.E. Raichev and F.T. Vaskcfl 

Institute of Semiconductor Physics, 
National Academy of Sciences of Ukraine, 
Prospekt Nauki 45, Kiev, 03028, Ukraine 
(Dated: February 6, 2008) 

Abstract 

Transient magnetotransport of two-dimensional electrons with partially-inverted distribution ex- 
cited by an ultrashort optical pulse is studied theoretically. The time-dependent photoconductivity 
is calculated for GaAs-based quantum wells by taking into account the relaxation of electron distri- 
bution caused by non-elastic electron-phonon interaction and the retardation of the response due 
to momentum relaxation and due to a finite capacitance of the sample. We predict large-amplitude 
transient oscillations of the current density and Hall field (Hall oscillations) with frequencies cor- 
responding to magnetoplasmon range, which are initiated by the instability owing to the absolute 
negative conductivity effect. 
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I. INTRODUCTION 



The transient process under an abrupt turn on of the electric current through a conduct- 
ing sample is described by a simple exponential dependence if the applied voltage is fixed, i.e. 
the load resistance is small and the circuit effects are not essential. The characteristic tempo- 
ral scale of this process is determined by the momentum relaxation time, which depends on 
the average energy of electrons and on the mechanisms of electron scattering. 1 ' 2 Such kind of 
exponential relaxation of electric current in pure bulk Ge has been demonstrated more than 
30 years ago. 3 To investigate the transient processes, one can use ultrafast photoexcitation 
of carriers instead of an abrupt switching of the applied voltage. Since this excitation cre- 
ates non-equilibrium electrons inside the conduction band, the temporal dependence of the 
current contains a slow component reflecting the energy relaxation of these electrons owing 
to quasielastic scattering by acoustic phonons. This slow energy relaxation, corresponding 
to the temporal scale which is much larger than the momentum relaxation time, takes place 
for the electrons in the passive region (i.e. for the electrons whose energies e are smaller 
than the optical phonon energy hu ), because the interaction of these electrons with optical 
phonons can be neglected at small temperatures T <C Tiu . The slow temporal dependence 
of the transient photocurrent in this case occurs for the same reason as the dependence of 
the photocurrent on the excitation energy under stationary photoexcitation. 4 

If the energies of excited electrons are close to hu , see the initial distribution B in 
Fig. 1 (a), the absolute negative conductivity (ANC) should take place, because the non- 
equilibrium distribution f £ of electrons in the passive region becomes inverted (df e /de > 0) 
in a certain interval of energies near the upper boundary of this region. Though such a 
possibility has been discussed a long time ago 5 for the regime of stationary photoexcitation 
in bulk samples, the ANC effect has not been detected so far in this regime. The absence 
of the ANC under a stationary photoexcitation is described by accumulation of low-energy 
electrons with time (owing to slowness of interband recombination) so that the relative 
contribution of the inverted high-energy part of electron distribution to the conductivity 
becomes non-essential. In addition, the Coulomb scattering of high-energy electrons by the 
low-energy ones leads to a rapid broadening of the initial narrow distribution of photoexcited 
electrons in the energy space, thereby decreasing the contribution of the inverted part of 
electron distribution. To date, the ANC has been realized by means of electron heating by 
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electric field in many- valley semiconductors 6 owing to intervalley redistribution of electrons, 
or when acting by microwave radiation on two-dimensional (2D) electrons in a quantizing 
magnetic field. 7 Recently, see Refs. 8, 9 and Ch. 11 in Ref. 10, it has been shown that the 
transient ANC, which exists during a time interval of the order of momentum relaxation 
time, can be achieved by ultrafast photoexcitation of electrons with energies close to hu . 
The theoretical description of this effect 9 has been based on the kinetic theory taking into 
account temporal non-locality of the response on the scale of momentum relaxation time and 
involving a detailed consideration of inelastic scattering of high-energy electrons by acoustic 
phonons. In this paper, we study the influence of classical magnetic fields on the transient 
ANC in 2D samples with the geometry of a long Hall bar, see Fig. 1 (b). 

The main feature of the transient magnetotransport under consideration is the appearance 
of temporal oscillations of the longitudinal conductivity and transverse electric field (Hall 
field), whose frequencies are of the order of the cyclotron frequency u c (Hall oscillations). 
Weak oscillations of this kind should be always present in the transient response because of 
the retardation of charge accumulation on the sides of the Hall bar. The existence of the 
transient ANC leads to the instability which dramatically modifies the transient oscillations. 
In the initial moments of time, when the longitudinal current flows in the direction opposite 
to the applied field, the Hall field increases in the direction opposite to its equilibrium one, 
because the sign of the Lorentz force is changed in the ANC regime. For the same reason, this 
increase is exponential: the charge accumulation on the sides leads to further enhancement 
of this accumulation. In the subsequent moments of time, when the partial inversion of the 
electron distribution is no longer sufficient to provide the ANC, the system starts decharging, 
and, since the system posesses a certain inertia, there appear large-amplitude oscillations 
of the Hall field as well as of the longitudinal current (this current is coupled to the Hall 
field). Both the current and the Hall field can change their signs in the process of the 
oscillations. The damping of such large-amplitude oscillations proceeds slower, which makes 
them persistent on a nanosecond time scale. 

The paper is organized as follows. In Sec. II we derive general equations for transient 
magnetotransport of electrons, which describe temporal dependence of the longitudinal cur- 
rent and Hall field. Section III is devoted to a simple model which makes it possible to solve 
these equations analytically and to describe the main features of the transient response. In 
Sec. IV we present the results of numerical calculations involving a detailed consideration of 
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FIG. 1: (a) Initial electron energy distribution /i for the cases of excitations away from the 
optical phonon energy (A) and near the optical phonon energy (B). (b) Hall bar geometry and 
electric currents in the presence of a magnetic field H directed perpendicular to the 2D plane. 

the evolution of electron distribution. The discussion of the assumptions used and conclud- 
ing remarks are given in the last section. The Appendix A contains the expressions for the 
transition probability and relaxation rate of 2D electrons interacting with acoustic phonons. 
The Appendix B contains the formalism describing the retardation of charge accumulation 
at the sides of the Hall bar. 



II. TEMPORAL RESPONSE 

We describe the response of photoexcited electrons to an electric field E t by representing 
the distribution function f pt , which depends on the electron momentum p and time t, in 
the form f pt = f £t + Af pt , where f £t is the symmetric part of this function and Af pt is the 
antisymmetric contribution induced by the field. Under the approximation of weak electric 
field, when the heating of electrons by the field is neglected, the symmetric part, which 
describes the energy distribution of non-degenerate electrons, is governed by the kinetic 
equation 

^ = Jac{ft\e) + J opt {f t \e). (1) 
The collision integral due to acoustic-phonon scattering, J ac , can be written as 

Jac (ft I e) = / de 1 [W £ , £ f e , t - W ££ ,f £t \ , (2) 
Jo 

where p S2 D = rn/rch 2 is the density of states of 2D electrons with the effective mass m. The 
scattering probabilities W £ > £ and W ££ < are determined by the deformation-potential (DA) and 
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piezoelectric (PA) interactions of electrons with acoustic phonons, see Appendix A. These 
probabilities satisfy the requirement of detailed balance, W e > e = W ee >exp[(e' — e)/T], where 
T is the phonon temperature. The collision integral J opt , which describes the interaction of 
non-degenerate electrons with dispersionless optical phonons at T <C Tiuj 0) can be represented 
in a similar form: 



where v is the rate of spontaneous emission of optical phonons by 2D electrons (see, for 
example, Ref. 10). 

Since the rate v Q is typically much larger than the rate of acoustic-phonon scattering, any 
electron appearing in the active region (e > hoj ) after photoexcitation or after acoustic- 
phonon absorption instantaneously relaxes to a state in the passive region, with the energy 
e — khu , where k is the number of emitted optical phonons. In this approximation, the 
kinetic equation (1) can be considered for the passive region only. To carry out such a 
transformation, we first rewrite the term describing departure of electrons from the state 
e in Eq. (2) as p^D Jo"" de' J2kLo W ££ > + khu> fet- Since the active region is empty owing to 
rapid emission of optical phonons (f et = at e > hcu ), this term is considered in the region 
e < Tiuj only. For the same reason, the term corresponding to arrival of electrons at the state 
e in Eq. (2) is written as p^D J^" de'W £ < £ f £ > t . This term describes transitions of electrons 
both to the states with e < hu Q and to the states with e > Tiu . As explained above, in the 
latter case the electrons instantaneously relax to the states with the energies e — khu in the 
passive region. Therefore, in the presence of rapid spontaneous emission of optical phonons 
the arrival term takes the form p^D Jq W ° de' Y^kLo W £ ' £+ khw fe't, where e < hu . Finally, 
since W ££ > becomes exponentially small at e' — e > hu , see Appendix A, one should retain 
only the terms with k = and k — 1 both in the arrival and in the departure terms. As a 
result, Eq. (1) is reduced to the following form: 



The term with W £ i £+ y lu>o in this equation becomes essential only if e' is close to Hui and e is 
close to zero. Similarly, the term with W ££ u r r llJ)o becomes essential if e is close to Tiuj and e' 
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(3) 




(W ££ , + W ££/+huJo )f £t ] . 



(4) 



is close to zero. 

Equation (4) should be solved with the initial condition f £t=0 = f^ ex \ where is 
determined by the excitation pulse. If the initial ultrafast excitation creates electrons with 
the distribution Fjf x \ we have f( ex ^ = 6(hi0 o — e) Y^kLo Fe+khu > wnere the terms with k ^ 
describe the electrons instantaneously transferred to the passive region via optical phonon 
emission. Note that Eq. (4) satisfies the density conservation requirement implying that the 
electron density pe,2D Jo"" def £t does not depend on time and equal to the excited density 

n ex = p,2D^°d^ ex) - 

The weak antisymmetric part Af pt , in the presence of a stationary magnetic field H 

(E_LH), is governed by the linearized kinetic equation 

where e is the electron charge, c is the velocity of light, and v = p/m is the electron velocity. 
The momentum relaxation rate on the right-hand side of Eq. (5) is the sum of the rate of 
quasielastic scattering of electrons by acoustic phonons, iA ac ), (see Appendix A) and the 
rate of spontaneous emission of optical phonons, v Q[e — Tiu ). The exact solution of Eq. 
(5) is given by 

A/ pi = e jf dt' e-*Wv • K* (-^j , (6) 

Kw = E t , coscj c (t - if) + ^ c X Et '^ sin w c (t - t'), 

u c 

where uj c = \e\H/mc is the cyclotron frequency vector. 

The current density is given by the standard formula, ji = (2/L 2 ) Z) p vA/ pt , where L 2 is 
the normalization square. Using Eq. (6) and performing the averaging over the angle of p, 
we write j t as 

it = ?£mf dt , Ka ,r deee -,t<-vL?M. (7) 

m Jo Jo \ as J 

The linear response of electron system to the electric field E t is described by the time- 
dependent conductivity tensor a tt > introduced according to the non-local relation j t = 
Jq dt'a U '~E t > . The diagonal and non-diagonal components of this tensor, of J and of/, are 

coso; c (t — t') 

(8) 

sina; c (t — t') 
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The contribution of the active region is neglected in this equation, because this region is 
depleted of electrons owing to rapid emission of optical phonons. Accordingly, the scattering 
rate v e standing in Eq. (8) is equal to the acoustic-phonon scattering rate iA ac ) calculated 
in Appendix A. 

Below we consider a sample of Hall bar geometry, a 2D strip of width d with the in-plane 
current density ] t = (jf ", jf- 1 ), where || and _L components of } t are referred to the coordi- 
nate system associated with the geometry of the Hall bar (Fig. 1 b). The transverse current 
density jf 1 - is not equal to zero and describes the transient process of charge accumulation 
on the sides (edges) of the Hall bar. Owing to near-edge localization of the magnetoinduced 
charges and current continuity, one can use the homogeneous current vector. Under the 
assumption of high resistance of the photoexcited electron gas in comparison to the load 
resistance of the circuit, we have E t = (E 6 ^, Ef 1 -), where the longitudinal field i? 6 " is deter- 
mined by the applied voltage and remains time-independent. The Hall field Ef 1 - depends 
on time because of the charge accumulation process mentioned above. The components of 
the current density vector are written through the components of conductivity tensor (8) as 
follows: 



where we have introduced the effective conductivity a t . To obtain a closed equation for 
the Hall field, one should describe the charge accumulation at the sides of the Hall bar. This 
leads to the approximate equation 



where Cg± = aed is the effective capacitance proportional to the dielectric permittivity e, 
and a is a small numerical factor. The derivation of Eq. (11) and the estimate of a are 
given in Appendix B. 

In summary, to describe the linear response of the system, one has first to solve Eq. 
(4) and determine the energy distribution f £t . Then, of J and of/ are calculated accord- 
ing to Eq. (8). Using them in Eq. (9) and applying Eq. (11), one finds the Hall 
field Ef-*-, which is proportional to the time-independent longitudinal field E B K Finally, 




(10) 



(9) 
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the longitudinal current is expressed through the effective conductivity from Eq. (10): 

ILL ANALYTICAL APPROACH 

Before presenting the results of numerical solution of Eq. (9), we discuss an approximation 
which allows one to understand essential features of the time-dependent response by means 
of analytical consideration. First of all, we neglect the energy dependence of the momentum 
relaxation time, replacing v e in Eq. (8) by a constant v = Uf luJo . Note that the calculated 
energy dependence v e in the interval < e < hu Q is not strong, except for the low-energy 
region (see Fig. 2). The integral over energy in Eq. (8) in this case is taken by parts, with 
the result 
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FIG. 2: Energy dependence of momentum relaxation rates for GaAs quantum well of width 10 
nm for T = 4.2 K and T = 20 K. The dashed and dotted lines show the partial contributions of 
DA and PA scattering mechanisms, respectively. 
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The evolution of the energy distribution enters the conductivity tensor (12) through the 
dimensionless function g t , which depends on the distribution function at the boundary of 



the passive region. After substituting the expression (11) to Eq. (9) with a J, and af t f given 
by Eq. (12), it is convenient to differentiate the equation obtained over t twice. As a result, 
Eq. (9) is reduced to the differential equation 

d 3 E^ n d 2 E? ± . 2 2 ^ 2 JE^ 

with the initial conditions £ t B± = dE? ± /dt = d 2 E? ± /dt 2 = at t = 0. We have introduced 
a characteristic frequency 

= /Z^I (14) 

determined by the capacitance and electron density. The longitudinal current, given in the 
integral form by Eq. (10), can be expressed through the derivatives of E^ with the use of 
Eqs. (9), (10), (11), and (12): 

7 B " - e2n ™ + u*?t + tfa eA (15) 

Jt ~ mWu e { dt 2 +U dt +ngtE * J' (15) 

If the electrons are excited considerably below the optical phonon energy, see the initial 
distribution A in Fig. 1 (a), one has g t = 1. In this case Eq. (13) is solved analytically. The 
solution shows a three-mode behavior according to 

E f± = (l _ Cie sit _ C2e s 2 t _ C3e s 3 t^ 5 (16) 

S 2 S 3 S1S3 
c l = 7 77 7 > c 2 



(si - s 2 )(si - s 3 ) ' (s 2 - si)(s 2 - s 3 ) 

_ S1S2 

(s 3 - si)(s 3 - s 2 ) 

where Si_ 3 are the roots of the characteristic equation s 3 + 2us 2 + (uo 2 c + u 2 + ft 2 )s + z/ft 2 
Under the approximation v 2 <C lj 2 + ft 2 , the solution (16) is rewritten as 



V 



l-exp(-g^) (17) 



ft 2 z/ ( co 2 + n 2 /2 



ex P !L_ n2 ^ sin(ft c t) 



ft? \ ft? 



where ft r = ^^Tft 2 . The expression (17) describes the increase of the Hall field from to 
its equilibrium value — E^uj c /u with a characteristic time (ft c /ft) 2 z/ _1 and weak oscillations 
of this field with the frequency ft c . The oscillations are exponentially damped on the time 
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scale of z/ _1 , though the damping is suppressed at (oo c /VL) 2 <C I. In the case of oo c 3> 0,, 
which still can be realized in the classical region of magnetic fields if the excitation density 
is low enough, the increase of the Hall field appears to be much slower than the damping of 
the oscillations. The longitudinal current shows a similar evolution, which is obvious from 
the relation (15). 

If the electrons are excited close to the optical phonon energy (see the distribution B in 
Fig. 1 (a)), the function g t is not equal to unity and can be negative at the initial moments 
of time, owing to the partial inversion of the electron distribution. As the excited electrons 
relax and go away from the boundary of the passive region, g t changes its sign from negative 
to positive at some instant t = t and approaches 1 at t — > oo. Although Eq. (13) cannot be 
solved analytically in the general case, the basic features of the response can be determined 
by using the model step-like function 

go, t<t 

9t= \ , (18) 

1, t > t 

where go is a negative constant. Substituting the expression (18) into Eq. (13), one can find 
that at t > t the solution (16) is valid again. However, the coefficients Ci_ 3 should be found 
by means of matching this solution to the solution at t < to, which has the form 



= _e*\\ ^ (l - d ie pit - d 2 e P2t - d 3 e P3t ) , (19) 



= - E"^' 

t<t 

and pi_ 3 are the roots of the equation p 3 + 2up 2 + (uo 2 c + v 2 + £l 2 go)p + v& 2 go = 0. The 
coefficients c?i_3 are expressed through pis in the same way as the coefficients Ci_3 are 
expressed through Si_3, see Eq. (16). The rules of the matching are derived from integration 
of Eq. (13) across the point t = t and imply continuity of the Hall field Ef 1 - and its first time 
derivative, while the second derivative has a finite step, d 2 Ef ± /dt 2 t _° + Q = —(1 — 9o)Ef^ tQ , 
which provides the continuity of the current given by Eq. (15). It is essential that at least 
one of the rates pis has a positive real part, which describes exponential increase of the 
Hall field in the interval t < to- This is a manifestation of the instability generated by the 
ANC effect. It is also important that the sign of the increasing Hall field is opposite to its 
equilibrium sign because of inversion of the direction of current in the ANC interval. The 
strong enhancement of the Hall field at t < t initiates large-amplitude oscillations of this 
field (Hall oscillations) in the region t > t . The oscillations of the longitudinal current are 
also dramatically enhanced. 
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FIG. 3: Evolution of the effective conductivity o\ e ^ and Hall field Ef 1 - calculated within the 
approximation described by Eqs. (12) and (18) at Q,/u = 5, to = v~ l , and go = —1 for two 
values of cyclotron frequency, uj c /u = 2 and uj c /v = 6. The dashed lines show the corresponding 
evolution for the case g t = 1, when the excited electrons appear considerably below the optical 
phonon energy. 

Figure 3 demonstrates some examples of temporal dependence of Ef 1 - and of 6 ^ (ex- 
pressed in units of a = e 2 n ex /mv) calculated within the analytival approach described 
here. We have chosen Q/u = 5, t = u^ 1 and g = —1. For comparison, we also plot the cor- 
responding temporal dependence at g t — 1, when the electrons are excited below the optical 
phonon energy. The oscillations in this case also exist, but they are weak and superimposed 
on a smooth exponential relaxation dependence. On the other hand, the oscillations gen- 
erated by the ANC instability are strong and remain visible even at vt ~ 10. Estimating 
v ~ 12 ns™ 1 , see Fig. 2 for 4.2 K, one can conclude that these oscillations persist over a 
nanosecond interval of time after the excitation. The frequency of the oscillations increases 
with the increase of the magnetic field and is estimated as Vl c = Joo^ + VL 2 . According to 
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Eq. (14) and to the estimate of Ce±, the frequency is of the order of 2D plasmon fre- 
quency uj q at wavenumbers q corresponding to the inverse width of the Hall bar, q ~ 1/d. 
For this reason, Q c is identified with a 2D magnetoplasmon frequency. The appearance of 
the large-amplitude Hall oscillations can therefore be considered as excitation of 2D magne- 
toplasmons owing to the ANC instability. The amplitude of the oscillations exponentially 
increases with the increase of the frequency Q and with the increase of the absolute value 
of g . We also note that the parameters used in the calculations are physically reasonable. 
Indeed, if v ~ 12 ns" 1 (Fig. 2), then Vt/v = 5 corresponds, for example, to C 6J _ ~ 0.1 cm 
and n ex ~ 10 11 cm~ 2 (or Cr± ~ 0.01 cm and n ex ~ 10 10 cm -2 ) for GaAs wells, while the 
ratios oo c /v = 2 and oo c /v = 6 approximately correspond to the classical magnetic fields of 
0.01 T and 0.03 T, respectively. 



IV. NUMERICAL RESULTS 



The model consideration given above ignores a detailed evolution of electron distribution 
after the photoexcitation. Below we present the results of a careful consideration based on 
the numerical solution of Eq. (9) with ofj and of/ given by Eq. (8). We have used standard 
material parameters of GaAs, which can be found, for example, in Refs. 1 and 11 (see also our 
paper, Ref. 9). To find the electron distribution f £t from Eq. (4), we assume that the optical 
pulse creates electrons with a Gaussian energy distribution oc exp[— (e — e ex ) 2 /A 2 ] of 
a characteristic half- width A, centered at the excitation energy e ex . Assuming, for example, 
that e ex = Tiuj and (A/hu ) 2 <C 1, we obtain the initial electron energy distribution in the 
passive region in the form of two half-peaks also shown in Fig. 1 (a): 

oc exp[-(e - hou ) 2 /A 2 ] + exp[- £ 2 /A 2 ], (20) 

see the discussion after Eq. (4). The numerical solution of Eq. (4) has been carried out by 
iterations in the time domain. The evolution of the electron distribution, calculated for 10 
nm wide GaAs quantum well at 4.2 K, is shown in Fig. 4. This evolution is similar to that 
calculated for bulk samples in Ref. 9. 

Figure 5 shows the temporal dependence of the effective conductivity crf e ^ and Hall field 
Ef 1 - for 10 nm wide GaAs quantum well at 4.2 K, calculated for the same parameters of 
the excitation. The effective conductivity is expressed in units of <r defined in Sec. III. The 
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FIG. 4: Temporal evolution of electron energy distribution for the case of initial distribution (20) 
with A = 0.1 hco , calculated for 10 nm wide GaAs quantum well at 4.2 K. The distribution 
functions are plotted with the time interval of 0.1 ns and normalized by M = n ex [p^2D^o] 1 ■ 

temporal dependence of <j\ e ^ at zero magnetic field is shown by a dashed line, and is similar 
to that calculated for bulk samples in Ref. 9. The evolution of a\ e ^ and Ef 1 - appears to 
be very sensitive to the cyclotron frequency and characteristic frequency Q [see Eq. (14)] 
because of the initial exponential increase of the current and Hall field. We have chosen these 
frequencies in such a way that the absolute values of of 6 ^ and Ef 1 - in Fig. 5 are not too 
large. By varying the parameters, one can obtain a very strong (several orders of magnitude) 
enhancement of <rf e ^ and Ef L , but the qualitative picture of the damped oscillations remains 
the same. The increase of the magnetic field leads to increasing frequency of the oscillations, 
while the amplitude of the oscillations decreases and the relaxation slows down. The increase 
of the excitation density, which leads to increasing frequency Q, exponentially increases the 
amplitude of the oscillations. 12 This influence of the parameters on the evolution is also 
described by the simple model investigated in the previous section. 

The general features of the evolution are not modified if the electrons are excited by a 
shorter optical pulse, which results in energy broadening of the initial distribution. Figure 6 
shows the temporal dependence of af e ^ and Ef 1 - calculated for the case of initial distribution 
(20) with A = 0.2 Tiuo . A comparison of this figure to Fig. 5 also demonstrates the increase 
of the oscillation frequency and a suppression of the damping as a result of the increased Q 
(see the discussion of Eq. (17)). 
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FIG. 5: Evolution of the effective conductivity crf e// and Hall field at Q = 50 ns 1 in the 
magnetic fields H = 0.015 T and 0.02 T, for the case of initial distribution (20) with A = 0.1 Uuj . 



The dashed line shows of e ^ (multiplied by 10) for zero magnetic field. 



V. CONCLUSIONS 



We have described the transient classical magnetotransport of electrons in a long Hall 
bar after ultrafast interband photoexcitation, and calculated the temporal dependence of the 
current and Hall field. Investigating the modification of the response due to accumulation 
of the edge charges forming the transverse (Hall) voltage, we have found the oscillations 
of both the current and the Hall field. The amplitude and duration of the oscillations 
are dramatically enhanced if the energies of the excited electrons are in the vicinity of the 
optical phonon energy. This is the case when the oscillations are triggered by the instability 
associated with the partial inversion of electron distribution (the ANC effect). Although the 
numerical calculations have been carried out here for non-doped GaAs quantum wells, similar 
effects should be expected for non-doped bulk samples, because the qualitative features of 
the energy relaxation and non-local temporal response of electron system do not depend on 
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FIG. 6: Evolution of the effective conductivity of e// and Hall field Ef 1 - at n = 80 ns" 1 in the 
magnetic fields H = 0.015 T and 0.03 T, for the case of initial distribution (20) with A = 0.2 Uuj . 
The dashed line shows a t e for zero magnetic field. 

dimensionality. 

Now we discuss the assumptions used. The main approximation is the consideration of 
electron scattering by phonon modes only. The elastic scattering by inhomogeneities can 
be neglected in the case of non-doped quantum wells with high-quality interfaces. In any 
case, it is not difficult to include this scattering into consideration, because it does not 
contribute to the energy relaxation (see Eq. (4)), and can lead only to an increase in the 
momentum relaxation rate standing in Eq. (8). However, the damping of the oscillations in 
this case becomes stronger, and it is always better to avoid this additional scattering by using 
clean samples. A more important restriction is the neglect of electron-electron interaction, 
which should dominate the energy relaxation of photoexcited 2D electrons at the densities 
n ex > 10 10 cm" 2 . Since the electron-electron interaction leads to a faster relaxation, it is 
expected to shorten the time interval where the exponential increase of the current and Hall 
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field takes place. However, this interaction cannot cancel the large-amplitude oscillations 
of the current and Hall field. In this connection, we note that these oscillations have been 
examined in Sec. Ill on the basis of a model that is not sensitive to the detailed description 
of the energy relaxation. The only essential point is the existence of partial inversion of 
electron distribution during a finite interval of time after the photoexcitation, sufficient for 
realization of the ANC regime. For this reason, the position and energy broadening of the 
initial electron distribution, which are determined by the parameters of photoexcitation 
pulse, appear to be more important than the actual energy relaxation mechanisms. 

Let us discuss the other approximations. The general formalism has been based on 
the classical Boltzmann equation. This is sufficient for the subject of our study, because 
the intervals of times under consideration considerably exceed the quantum broadening 
times h/e, and the magnetic field is weak enough to neglect the Landau quantization. We 
have ignored the existence of holes created in the valence band by the optical pulse. This 
is possible because of smallness of the hole mobility, so the contribution of the holes to 
the transport can be neglected. To describe the momentum relaxation by electron-phonon 
scattering, we have used the elastic approximation. In quantum wells, this is suitable for 
description of electrons whose energies e are much larger than the characteristic energy hsu/a 
associated with phonon momentum normal to the layer (here s is the sound velocity and a is 
the well width). For typical parameters s ~ 5 x 10 5 cm/s and a ~ 10 nm, the energy hsir/a 
is around 1 meV. Therefore, the assumed condition is satisfied in our calculations for the 
photogenerated high-energy electrons which give the main contribution to the conductivity. 
Next, since we have neglected the transverse inhomogeneity of the currents and fields, the 
equations (9)-(ll) are not applicable for description of electrons in the vicinity of the Hall 
bar edges. Nevertheless, these equations are valid in the main part of the Hall bar, where 
the inhomogeneous corrections to the currents and fields are small (see Appendix B). 

We also note that in the ANC regime one should consider a possibility for the development 
of a spatial instability 13 both along and across the Hall bar. Uncovering the conditions for 
existence and properties of such an instability requires a special investigation based on 
the formalism of non-homogeneous kinetic equation. The spatially-inhomogeneous electron 
distribution owing to the ANC effect is essential in the stationary regime, and it is realized 
for 2D electrons under microwave excitation in the quantizing magnetic field. 7 Nevertheless, 
since we consider the transient response, one may expect that the results of this paper will 
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remain valid for the samples of small size, where the spatially-inhomogeneous distribution 
is not developed during the short interval of time corresponding to the exponential increase 
of the current and Hall field. 

Finally, we would like to point out that the Hall oscillations studied in this paper can 
be viewed as 2D magnetoplasmons with small wavenumbers determined by the width of 
the Hall bar. To detect them in experiment, it is necessary to have sub-nanosecond tem- 
poral resolution, which is attainable for standard all-electrical measurements. Owing to 
the strong amplification of the oscillations by the ANC instability, the excitation technique 
based on the ultrafast optical pump can be applied along with the conventional techniques 14 
of magnetoplasmon excitation in 2D layers. 



APPENDIX A: RELAXATION RATES 

Below we present the expressions for z/^( ac ) and W ee < determining the momentum and 
energy relaxation of 2D electrons under acoustic-phonon scattering. Both these quantities 
can be written as the sums of partial contributions caused by the DA and PA mechanisms 
of interaction. The relaxation rate z/?( ac ) is calculated in the elastic approximation: 



xS(e p -e p/ ) £ E |CaqT(2^aq + 1), (Al) 

i=RDA,PA X=l,t 

where e = e p = p 2 /2m, Q = (q, q R± ) is the phonon wave vector, q = (p — p')/h is the 
in-plane component of this vector, ip is the angle between p and p', N X q are the Planck 
occupation numbers of the longitudinal (A = I) and transverse (A = t) acoustic phonon 
modes, and |C*q | 2 are the matrix elements of interaction, averaged over the angle of electron 
momentum p in the 2D plane. This averaging is essential for the PA interaction, which is 
sensitive to orientation of the vector Q with respect to crystallographic axes. 15 In the explicit 
form, 

^Ty-T^-^ (A2) 



and 



2ps x 



lr &(PA) {2 _ h(eh 14 ) 2 . ( . s 
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where s\ are the sound velocities, p is the material density, V is the deformation-potential 
constant, and hu is the piezoelectric coefficient. The polarization factors A\q for (100)- 
grown 2D layers are 16 

9 ql , g 4 , A qLq 2 1 q 6 , . „, 

A «* = YW' A ^ = i ~W + 2W (A4) 

Finally, M(q^±) = |(0|e* 9±z |0)| 2 is the squared matrix element of a plane-wave factor. This 
matrix element characterizes the interaction of 2D electrons with 3D phonons and depends 
on the confinement potential determining the ground state of 2D electrons, |0). We apply the 
model of a deep square well of width a leading to the expression 2 M(q R± ) = [sm(x)/x] 2 /[I — 
(x/ir) 2 ] 2 , where x = q&±a/2. 

The introduction of the transition probability W e£ > entering Eqs. (2) and (4) implies 
averaging of the kinetic equation over the angle of momentum p, according to df £ t/dt = 
(2n)~ 1 J 27r d(f p E P '(Wp'p/e't — Wpp'/et)) where W pp > is the probability of electron scattering 
from the state with momentum p to the state with momentum p'. As a result, 

1 sgn(e — e') 



W FF , 



x 

i=DA 



~~ h 2 l-exp[(e> -e)/T] 



e e/;>j^, (A5) 

DA, PA \=l,t 71 S A?B± 



where |C^q | 2 and M(g 6J _) are already defined above. However, the wave numbers Q, q, 
and qe± standing in the corresponding equations should be now considered as functions 
of energies, phonon polarization, and scattering angle ip, according to Q — \e — e'\/hs\, 
q = \j2m(e + e' — cos<p) /%, and q & ± = \JQ 2 — q 2 . The integral over the scattering 

angle in Eq. (A5) must be calculated numerically. The upper limit of this integration, ip\, 
is determined from the requirement Q 2 > q 2 , which means (f\ = n at (y 7 ^ — \fe') 2 > 2ms\ 
and (fx = arccos (je + e' — (e — e') 2 /2ms 2 ]/2\/ee'^ at (y/e — \fe') 2 < 2ms 2 . In the case of 
very small energies, (y/e + Ve 7 ) 2 < 2ms 2 , the limit (p x should be set at zero and the integral 
gives no contribution. 



APPENDIX B: EDGE CHARGE ACCUMULATION 

This appendix contains a simple formalism describing inhomogeneous distributions of 
charges, currents, and fields in a long Hall bar for a time-dependent problem. The homo- 
geneous approximation, which is valid in the central part of the bar and leads to Eq. (11) 
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for the transverse current describing edge charge accumulation, follows from the general 
consideration. 

Consider a semiconductor structure containing a 2D Hall bar of width d and length 
L {—d/2 < y < d/2 and —L/2 < x < L/2) placed at a distance h below the surface 
of the medium with dielectric permittivity e. If the retardation effects are neglected, the 
electrostatic potential z) created by the 2D carriers satisfies the Poisson equation 



d d d 2 
—e(z)— + e(z)- 



. h ,-.,. h .,^. h/2 My,z) = -^e5n t (y)5(z + h), (Bl) 

where e(z) = e at z < 0, e(z) — 1 at z > 0, 5n t (y) = n t (y) — p is the excess density of 
free electrons, and p is the density of holes which is assumed to be homogeneous and time- 
independent (because of low mobility of holes, this is possible in the short time interval under 
consideration). The Poisson equation is two-dimensional because the bar is long, d <C L, 
and the electron density depends only on the transverse coordinate y. If e ^> 1, one can use 
the Newmann's boundary condition d$ t (y, z)/dz\ z= _ = at the surface and rewrite Eq. 
(1) in the integral form 

rd/2 
-d/2 

The other necessary equations are the continuity equation 



P rd/2 _ , 

$ t (y, *) = --/ dy> ln[(y - y'f + (z + hf] + ln[(y - y'f + (z - h) 2 } Sn t (y'). (B2) 
e J -d/2 1 J 



e^M + MiM = 0, (B3, 
at ay 

which relates the density to the transverse current, and the equation 

5n t (y) = ep mD [v t (y) - <p t (y)], (B4) 

where (pt(y) = ®t(y, ~h) is the electrostatic potential in the 2D plane and v t (y) is the 
electrochemical potential. Finally, there is a relation between the transverse current and 
electrochemical potential, 

j^(y) = fdt* (-aS^-Mv) + °u>E m ) , (B5) 



dy 

which is a generalization of Eq. (9): Ef, 1 - replaced by —dv t '(y)/dy in order to take into 
account the spatial inhomogeneity. 

Using Eqs. (B2), (B3), and (B4), one can exclude the electrostatic potential and electron 
density, and write 

- ~^y- + — r- L/2 dy K ^ y ] ^r~ = e p ^ D —> (B6) 
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where K(y,y') = In \y - y'\ + In ^J(y-y') 2 + (2h) 2 . Equations (B5) and (B6), with the 
boundary condition jf ± (±d/2) = 0, give a complete description of the spatial distribu- 
tion of transverse current and electrochemical potential in the long Hall bar. The spatial 
inhomogeneity is essential near the edges ±d/2, this follows from the fact that the trans- 
verse current goes to zero at y = ±d/2. On the other hand, in the main part of the bar 
one should expect nearly homogeneous currents and fields. If the bar width d is large 
in comparison to the Bohr radius h 2 e/e 2 m, the main contribution to the right-hand side 
of Eq. (B6) comes from the integral term, which is written after integration by parts as 
(2e 2 pB2D / e ) I-d/2^y'[^^(y->y')/^y\jt ± (y')- Let us search for the solution of Eq. (B6) in the 
form jt L {y) = jf ± + 5jt(y), where jf 1 - is the homogeneous part of the current and Sjt(y) 
is the inhomogeneous correction. Neglecting first Sj t (y) (zero-order iteration), we obtain 
jf 1 - = e{2[K(y, —d/2) — K(y, d/2)]}~ 1 [dv t (y)/dt], which is rewritten, under the reasonable 
assumption d ^> Ah, as 

•e± = i 9vt(y) ™ 

Jt A\n\(d/2 + y)/(d/2-y)\ dt ' 1 ] 

Generally speaking, this equation can be satisfied for all y if the transverse field is inho- 
mogeneous. In other words, v t (y) = —Ef^y + 5v t (y), where 5v t (y) is the inhomogeneous 
correction to the electrochemical potential. However, in the central part of the sample, 
where In | (d/2 + y)/(d/2 — y)\ ~ Ay/d + 0[(y/d) 3 ], the inhomogeneous correction 5v t (y) 
can be neglected and the electrochemical potential is expressed through the homogeneous 
transverse field, v t (y) ~ —Ef ± y. Equation (B7) in these conditions is reduced to Eq. (11) 
with the effective capacitance Cqj_ = ed/16 (so the coefficient a introduced after Eq. (11) is 
equal to 1/16), while Eq. (B5) is reduced to Eq. (9). Using these solutions, one can write an 
equation for the inhomogeneous correction to the current (first-order iteration), which shows 
that the relative contribution of Sj t (y) is small, d~ x J^ 2 dy5j t (y) jf- 1 , this smallness is 
of numerical origin. Therefore, the homogeneous approximation described by Eqs. (9) and 
(11) is justified. 
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